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RESUMEN 

Presentamos un metodo que utiliza modelos de discos para determinar masas 
de estrellas T Tauri a partir de fotometri'a JHK. Encontramos que el exceso infrar- 
rojo producido por el disco se comporta de forma similar a un vector de extincion 
cuando son graficados en diagramas color-color y color-magnitud, pero son lineal- 
mente independientes. Mediante la manupulacion de estos vectores como si fueran 
una base de un espacio vectorial, Uevamos a cabo una transformacion de coorde- 
nadas que nos permitio encontrar la masa de la estrella central para una edad dada. 
Para probar este metodo basado en vectores principales (PV) , comparamos los val- 
ores de masa obtenidos (Mpy) en 14 objetos T Tauri de la region de Taurus- Auriga 
con los valores de masa calculados por metodos dinamicos. Tambien analizamos 4 
sistemas multiples (GG Tau, FO Tau, FS Tau y V733 Tau) y determinamos la masa 
de cada componente del sistema. Se analizo el efecto que se produce al utilizar difer- 
entes modelos evolutivos y la variacion de otros parametros en la determinacion de 
masas este metodo. 

ABSTRACT 

We present a method that uses disk models to determine stellar masses for 
Pre-Main sequence T Tauri stars from JHK photometry. We find that the infrared 
excess produced by the disk behaves in a similar way to the extinction vector 
when plotted on color-color and color-magnitude diagrams, but that it is linearly 
independent. Employing these vectors as a basis of a vector space we carry out 
a coordinate transformation that allows us to find the mass of a central star for 
a given age. To test this Principal Vectors (PV) method we compare the mass 
values (Mpv) obtained for 14 T Tauri objects in the Taurus- Auriga region with 
mass values known by dynamical methods. Further on, we analyzed 4 systems with 
multiple components (GG Tau, FO Tau, FS Tau, and V733 Tau) and determined 
the mass for each component of the system. We analyzed the effect of using different 
evolutionary models and other parameter values. 

Key Words: STAR FORMATION: GENERAL STAR FORMATION: 
T TAURIS, PHOTOMETRY, PRE-MAIN SEQUENCE 
MASS. 
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1. INTRODUCTION 

The only direct method to determine stellar 
masses is through the analysis of dynamical param- 
eters. This can be accomplished by measuring peri- 
ods of orbiting components, radial velocities of un- 
resolved companions, or from the kinematics of a 
circumstellar disk. Obtaining the necessary data is 
often impossible, or else difficult, and so dynamical 
methods have been used only for a few pre-main se- 
quence (PMS) stellar objects. As a consequence of 
this, most of the mass values are obtained by sec- 
ondary methods. 

The usual secondary method is to place the ob- 
jects on a Temperature-Lu minosity diagram alon g 
with evolutionary models (jCohen fc Kuhil Il979l ). 
Many circumstances encumber this task. First, tem- 
perature could be characterized from the colors of 
the objects. But colors are affected by circumstellar 
disks. The disks obscure optical observations and 
produce infrared (IR) excess, also the accretion onto 
the surface of the star increases the UV part of the 
spectrum and veils the optical, and a reddening is 
produced by extinction from the dusty medium. 

Therefore, temperature is better determined 
from the spectral typ e of the object, obtained 
throu g h spectroscopy (ILeyine et al. 1 120061: IdJ 
20051 : iLuhman et"al1 l2005t iLuhman et all l2003t 



Briceno et al. 20021: iLuhman [20001: ILuhman et al 



20001: iLuhmanI Il999l : iLuhman fc Riekd Il999l: 



Luhman et al.lll998f ). To determine the luminosity 



of the source, extinction corrections have to be ap- 
plied. These corrections can also be determined from 
the spectral type of the source by comparing the ob- 
served colors with the intrinsic ones. In doing this 
task, the use of colors where either the UV excess due 
to accretion or the IR excess due to disk emission are 
usually avoided because both effects contaminate the 
colors. Thus, co lors involving the bandpases R, I an d 



J are preferred (jLuhman et al 



2003t lLuhman|l2000l) 



Among several improvements. iHartigan et al.l (| 19911 ) 

have refined the method by modeling the veiling to 
correct the spectrum and so, by this way, he obtains 
a better estimate of the spectral type and the extinc- 
tion corrections. However, the acquisition of spectro- 
scopic data is a difficult task for faint objects in dusty 
environments, which constitute the vast majority of 
PMS stellar objects. This has lead many authors to 
explore the initial mass function (IMF) of clusters in 
the IR. 

The most used measure of the IMF is the 
K luminosity function (KLF) because it is likely 
to produce a fait hful representation o f the true 
IMF of a cluster (jMuench et al. 1 12002|) . This is 



so because the IR excess is only a small correc- 
tion when large KLF bins are employed. However, 
the K luminosities of objects in clusters are af- 
fected by extinction. Although this problem some- 
times is not accounted for, most authors either 
correct for ex tinction using the color-color dia- 
gram fC-C PI (iLeistra et al-lbood: iHerbig fc Dahm 
20061 : iMcCabe et al.l l2006t iFigueredo et all l2002t 



Bontemps et al. 2001; Luhman 2000l : Porras et all 
20001) or the color - magnitude diagrain (C-M D) 



(|Leistra et all 120061: IPorras et all l2000l ). On the 



other hand, the K luminosity is seldomly cor- 
rected by the IR e xcess emission (however see 
Figueredo et al.ll2002r ) and authors that worry about 
it, sometimes favor the use of th e JLF or HLF 
where the contamination is smaller dLiu et al.l 20031: 
Nielbock et 311120031: lB"randl et al.iri999l: lOasa et al 
1999f). 

But the IR excess emission has been clearly iden- 
tified and shown to h ave a definite ef fect i n C- 
C and C-M diag r ams. Lada fc Adamd (|l992h and 



Hillenbrand et al.l (|l992f ) show that objects of dif- 
ferent evolutionary st ates occupy definite regions of 
the C-C D (see also |Palla.lll993l) The classical T 
Tauri star locus ( Meyer et al.l 1997 ) is a well defined 
linear region in the C-C D. This locus is occupied 
by dereddened T Tauri stars from a Taurus sample. 
In their work, (jMever et al.lll997l ) displace T Tauri 
stars along the reddening vector by an amount deter- 
mined from the E(Rc-Ic) excess above the intrinsic 
colors known from the spectral type of each source. 
The resulting locus lies along a narrow line redwards 
of the main sequence colors corresponding to stars 
of spectral types around M. That is, the excess in 
the colors, in particular E(J-H) and E(H-K), are re- 
lated to each o ther and define an e xcess vector in the 
C-C D. Then, iMever et all (|l997l ) model the excess 
emission from accretion disks around K5 to M5 stars 
to show that this locus is indeed well explained and 
expected due to the existence of circumstellar disks. 
Therefore, it is proven that by adding an excess vec- 
tor and an interstellar extinction vector to a stellar 
photosphere, one can recover the position of a par- 
ticular T Tauri star on the C-C D. However, if the 
spectral type of the star is not known, the magni- 
tudes of the vectors can not be assessed. 

The existence of a vector like excess, this time 
for the K vs H-K diag r am, h as been found by 
Hillenbrand fc Carpenter ( 2000l ). They analyzed a 
similar sample of T Tauri stars in Taurus and 
found a linear relationship between the K excess 
emission and the excess in H-K. Again, the ex- 
cesses are determined from the dereddened colors 
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and magnitudes if the intrinsic colors are known 
from the spectral types of the so urces. In their work, 
Hillenbrand fc Carpentel (|200d ) e;o on to obtain the 
IMF of the trapezium cluster by using the vector 
sum of this excess vector and the interstellar extinc- 
tion vector. However, since they only used the C-M 
D, they had incomplete information to obtain the 
magnitudes of these two vectors, so they relied on a 
statistical approach that treats these unknown quan- 
tities as probability distributions. 

Therefore, it has been shown that IR excess emis- 
sion affects both near infrared (NIR) C-C and C-M 
diagrams in a similar way as the interstellar extinc- 
tion. However, the two diagrams have not been used 
simultaneously to obtain the magnitude values of IR 
excess and extinction. In this work we will explore 
this possibility. We will begin by analyzing the ef- 
fect of a circumstcllar disk on the NIR colors and 
magnitudes. Theoretical models to explain disks 
may turn out to be very com plicated and include 
many parameters (IJohns-KuU et al.i.2001, ). The first 
circumstellar disk models were assumed geometri- 
cally thin and optically thi c k (IShakura fc Sunyaev 
1973 ). Kenyon fc Hartmann ( 19871 ) realized that the 
thickness of the disk should increase radially, leading 
to a flared disk. This resulted ia a better fit to the 
spectral energy distribution (SED) for many T Tauri 
stars. This disk model assumed that the temperature 
was only dependent on the disk radius. The vertical 
temperature gr adients in th e disk were first taken 
into account bv lCalvet et al. I (ll99lLll99^. wit h fur- 



ther analysis by 'P'Ales sio et all (|l99'8ril999f) . and 
Chiang et all ([l997, ,199j). In these models the ab- 



sorption of stellar radiation occurs only in the top 
disk layers, so the heating of the inte rnal disk i s pro- 
duced by viscous dissipation. P'Alessio et al.l (|l998l . 
19991) calculated the spectral emission of a disk as- 
suming the sam e structure for t he gas and dust in the 
disk, whil e Chiang et al. ( 19991). ignored the gas con- 
tribution. iD'Alessio et al. I l|l998h show that if there 
is a low coupling between the gas and dust at the 
surface of the disk, the gas temperature decreases. 
In that case, the dust dominates the disk emission 
in classical T Tauri stars observations. 



Dullemond et al 



et al.l 
1997tl 



Chiang et all (1199 
up inner rim 



for 



(|200lD 
model to 
Ae/Be stars. 



extended the 
incl ude a puffed- 



Muzerolle et al. 



(|2003( ) showed that the emission of a puffed rim also 
fits the SED of T- Tauri stars in the 2-5/im region. 
The rim models were rather simplifie d treating them 



as a vertical walls ( as in models by ID'Alessio et al. 



20051 ). Rounded-o ff rims seem to be mor e consistent 
with observations ( Isella fc Natta I I2OO5I) . 



We will explore the effects of the disk excess emis- 
sion in C-C and C-M Diagrams using reprocessing 
models and more complicated accretion models. We 
will show that in certain cases the produced excess 
can be interpreted as a reddening vector. Then, we 
will show on the diagrams that these vectors point 
out in a different direction than the interstellar ex- 
tinction one, and so we will identify a pair of Prin- 
cipal Vectors (PV's) produced by the disk ( D) and 
interstellar extinction {X). These vectors can be used 
as a tool to determine stellar masses from JHK pho- 
tometry, without the need of spectral types. 

A set of PMS tracks at the cluster's age will 
be needed to read out the mass of the star. We 
will take advantage of rece nt work on evolution 
ary m odels 



like those of 
19971 ) (hereaft er DM97) 



D'Antona fc MazziteUi 



SOO) iBaraffe et al 



the Yale Group Y' 



. iPalla fc Stalher 

(1 99ai (hereafter PS99):ISiess et al.l (l2000|) (hereafter 
IT 



( 1998h (hereaft er BCAH98) and 



lYi et al.l (l200lh (hereafter F^qi). 

iHillenbrand fc Whitel ( 2004 ) reviewed these models 
listing the basic parameters employed by each one. 

In §2 we proceed with the calculation of the NIR 
excess emission of disks from reprocessing and ac- 
cretion models and obtain the PV for each diagram. 
Also, we will describe the procedure to obtain mass 
values with the PV method (Mpy) using synthetic 
evolutionary models for intermediate to low-mass 
stars. In §3, we test our method by applying it to 
a set of 14 objects of the Taurus-Auriga association 
(d^l4G pc) with previous mass determinations by 
dynamical methods. In addition, we will apply the 
method to multiple systems and determine the mass 
for each component. The conclusions are given in §4. 



2. PV METHOD 

We will begin by calculating the NIR effect of 
disks and show that the excess seen in CC and CM 
diagrams is similar to the behavior of the reddening 
vector produced by interstellar extinction. We will 
first illustrate this effect with a simple reprocessing 
disk via a simple and easy to manage calculation. 
The real situation, however, is more complex. Heat- 
ing of the disk is also accomplished via viscous dis- 
sipation during accretion, and a plethora of parame- 
ters contribute to the process. Thus, we present the 
NIR colors and magni tudes of accretion disk models 
(jD'Alessio et al.l 120051 ) in the following section. Ex- 
tinction is briefly described following that and then, 
we identify the principal vectors, describe how these 
vectors can be used to determine stellar masses and 
the expected limits of applicability of the method. 
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2.1. Excess emision vector for reprocessing disks 
models 

Passive disks models do not have an intrinsic lu- 
minosity and only re-emit the energy absorbed from 
the central star. Thus the superficial disk tempera- 
ture ca n be approximated by a power law expression 
(like in lNatta et"alll2000l (hereafter NOO)) hke 



Tiro) To {rD/R,y 



(1) 



where the rjj is disk radius, and i?, are the 
effective temperature and the stellar radius respec- 
tively. NOO find that the parameters Tq = 0.75T, and 
q=0.47 provide a good agreement with most observa- 
tions. We considered a collection of concentric rings, 
each one emitting as a blackbody in agreement with 
eq 1. The range of temperature goes from an inner 
radii with Tin^lSOO K (dust sublimation tempera- 
ture) to a maximum radii with Tout— WOK. Taking 
into account these considerations we calculate the 
predicted JHK disk flux. The total disk flux at a 
given frequency is given by. 



f[Dl= f B,(T{rD))dn, 
Jn 



(2) 



where B^{T{r£i)) is the Planck function and T{r£i) 
is the local temperature as a function of the disk 
radius. 

So, to carry out the integral getting the next ex- 
pression: 



2^-2/9 



-(2/<z) + l 
jl/^ - 1 



■dx. 



(3) 



where a — Airhv^ / qc^ , (3 = R^,/dobsj a-nd k — 
kT^/hv, the integral limits are given as a function 
of the disk temperature, being Tin and Tout the val- 
ues previously mentioned and 7 = k/hv. To obtain 
results from eq. 3, we only need to know the effective 
temperature of the central star T* = 7e//, the stellar 
radius i?*, and the frequency v. Once the disk flux is 
calculated, we add it to the stellar flux following the 
next expression that constitutes the star-disk system 
[s+d] 

f[s + d\, = f[MS*], + ef[D],, (4) 

where e represents several parameters such as the 
inclination angle to the line of sight and circumstellar 
dust properties. 

We use a set of intermediate to low mass 
main-sequence stars (Spectral Type from AO to 
M6) and construct the s+d system using effec- 
tive temperatures and bolometric magnitudes from 



Schmidt-Kalen (jl982f ) to calculate the stellar ra- 
dius and s tandard JHK photometric values from 
Bessell fc B rett (1988). In Figures la and lb we show 
the C-C (J-H vs. H-K) and C-M (K vs. J-K) dia- 
grams with the set of main sequence stars and the 
s-|-d systems . The results are simila r to those of oth- 
ers authors. Lada fc Adamd (1992), make use of C- 
C D to show how the IR colors are affected by the 
presence of circumstellar disks on intermediate to 
low mass stars, and achieved similar results. Another 
useful diagram is the Color-Magnitude. We use the 
K vs. J-K because with this combination we get the 
maximum contrast between the IR excess produced 
by the presence of the disk, which affects mostly K, 
and the interstellar extinction which has a greater 
effect in J. 

The main sequence is labeled with the spectral 
type of the star, the open circles represent the s-|-d 
model, the solid lines joining the two marks are the 
disk's vector direction and length, that is, the dis- 
placement produced by the presence of the disk. We 
show the case for e = 0.4 that corresponds to a 66° 
of inclination from the line of sight. From figures la 
and lb we perceive that it is possible to represent the 
IR excess produced by the disk as vectorial quantity 
in each diagram. This vector points in nearly the 
same direction for stars of all masses in the low to 
intermediate range. We only start to see a discrep- 
ancy in the direction for stars earlier than AO. For 
the rest of the spectral types, the median value of 
the displacements, represented by Z?, reproduces the 
mean effect of the disk. The length of the vector is 
however different for each case, but this only implies 
that an amount ||Z3|| of disk emission should be inter- 
preted differently for each spectral type. This same 
procedure was repeated for a collection of 10,000 ran- 
domly oriented viewing angles from to 90 deg . The 
resultant median vector D shows little variation in 
the direction. In The C-C D the slope of the excess 
emission vector is 0.43±0.107 and in the C-M D the 
slope is -1.097±0.02 . Thus, the obtained compo- 
nents of the excess vectors are: Z5cc=(0.731, 0.317) 
and Z)cm=(l-05, -1.15) in magnitude units. 

2.2. IR excess vectors for accretion models 

Following the modeling work of D'Alessio and 
collaborators (1998, 1999, 2001), they have pub- 
lished the SED's for a grid of 3000 accretion disks 
models that cover a wide range of parameters 



(jD'Alessio et al.l 120051 ). These models include accre- 
tion rates between 10~^ to 10~^ MQ/yr, onto cen- 
tral stars of 0.8 to 4 M© with temperatures from 
4000 to lOOOOK and ages of 1 to 10 Myr. Disk sizes 
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Fig. 1. (a)(J-H) vs. (H-K) and (b)(K) vs. (J-K) dia- 
grams used to determine the PV generated by the disk 
reprocessing vector {D) and interstellar extinction vec- 
tor {X). We use stars with spectral types from AO to 
M6, with a value of e = 0.4 that corresponds to a 66° of 
inclination from the line of sight . The directions of the 
vectors are shown on the top left part on each figure. 



range from 100 to 800 AU, and have inner radii that 
depend on the dust destruction temperature, taken 
around 1400 K, where the disk develops a wall whose 
scale height is also considered. Gas and dust are well 
mixed in the disk, with dust size distributions given 
by power laws with indexes of 2.5 and 3.5 and max- 
imum particle sizes of 1 /im to 10 cm. The systems 
are viewed at two inclination angles; cos(«) — 0.5 



and cos(i) = 0.86. 

From this collection of models, we have chosen 
those that provide a clear view of the star (cos(i) — 
0.86) so that the effects of extinction can be sepa- 
rated from those of IR excess. The effects of extinc- 
tion will be addressed in the next section. We also 
limit the following analysis to a particular age of 10 
Myr for clarity. Thus we are left with ~ 730 disk 
models. From the SED's of the star-|-disk systems we 
compute the magnitudes in the J. H and K ba nds 
(calibrating according to Besse ll fc Brett 1988h as 
well as the correspon ding magnitudes for ea ch cen- 
tral star, also given in lD'Alessio et al.l (j2005l ). 

In Figure 2 we present the C-C and C-M dia- 
grams. In each panel the colors and magnitudes of 
the central star are represented by circles labeled 
with the stellar masses (0.8, 1.1, 1.4, 1.6, 1.7, 1.9, 
2.0 and 2.3 Mq), and the crosses represent the cor- 
responding star-l-disk systems. 

The circles that corresponds to masses around 2 
Mq tend to pile up together, as the masses them- 
selves are very close together. However we have de- 
cided to keep the limits of the figure to allow a direct 
comparison with Figure 1. In the C-M diagram it is 
evident that the disk models for each star occupy a 
very definite region, lying along a unique direction 
to the upper right, that is, going redder in J-K and 
brighter in K. The median direction is indicated by 
an arrow, and it is very close to the direction found 
for simple reprocessing models (previous section) . In 
the C-C diagram, the disk models corresponding to 
each star also tend to occupy the upper right re- 
gion, that is, redder in both colors, although the di- 
rections show a little more dispersion than in the 
previous case (reprocessing models). In fact, the de- 
viations from a common direction of all vectors in 
this figure is 0.15 radian, while only 0.02 radian for 
the C-M diagram. Also note that the deviations are 
more important as the mass of the star increases, in 
the sense that the vectors curve upwards for shorter 
vectors. This fact will limit the range of applicability 
of the method that we are proposing, since we will 
rely on a common direction for this excess vector, 
as indicated by the arrow in the C-C figure. An- 
other important fact is that if we ratio the magni- 
tude of each excess disk vector in the C-M diagram 
to it's corresponding magnitude in the C-C diagram, 
we get a constant 1.73 ± 0.14, indicating that vec- 
tors behave accordingly in both diagrams. Thus, we 
may define the median disk vectors with components 
Z3ec=(0.767, 0.221) in the C-C D and D^n^ (1.014, 
-1.105) in the C-M D (actual median values have 
been arbitrarily multiplied by 13 to allow a better 
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Fig. 2. (a)(J-H) Vs(H-K) and (b)(K) vs. (J-K) diagrams 
used to determine the PV. The disk vector {D) is gener- 
ated by the median value of ~ 730 accretion disk models 
and the interstellar extinction vector is represented by 
X. 

comparison), which are indeed very similar to the 
ones derived for the simple reprocessing disk models 
of the previous section. 

2.3. Extinction Reddening Vector 

Normal interstellar extinction can be represented 
as a reddening vector with components [E{J — 
K),Ak] and [E{H - K), E{J - H)] in the C-M and 
C-C diagrams respectively. In Figures 1 and 2, we 
denoted by X the vector corresponding to Ay =1 
and normal colors given in iRieke fc Lebofskvl (Il985l) . 



It is our intention to use these vectors to quantify 
the reddening produced by the ISM, the circumstel- 
lar envelope and the disk itself, on the photosphere 
of the star and the NIR emitting region of the disk. 
However, we must proceed carefully in this task. It is 
known that R, the ratio of total to selective extinc- 
tion, attains larger than normal values in dense cores 
and dark clouds associated to star forming regions. 
This is thought to arise because the shielded condi- 
tions of high opacity regions favor the coalescence of 
small grains into larger particles. Larger particles in 
turn produce a more geometrical attenuation of vis- 
ible and NIR light, res ulting in a "g r ay" ex tinction. 

In a recent paper, iMoore et al. (|2005l) has ob- 
served the Hydrogen NIR recombination lines to- 
wards UC-HII regions, to show that the 2/im opac- 
ity becomes flatter than normal only when Ay > 25 
mag. If an appropriate value of R can be determined 
for these dense cores, the components of the redden- 
ing vector can be readily determined. 

The mid planes of circumstellar disks are also re- 
gions of high opacity that can affect the observations 
of a central star with abnormal extinction. However, 
D'Alessio et al.l ( 200ll) argue that this is only a prob- 
lem for a small (15%) fraction of T Tauri stars in 
Taurus (a typical star forming region). The increas- 
ing size in dust particles results in a growing critical 
viewing inclination at which the obscuring Ay > 30 
condition is reached. 

We can therefore assume that a normal extinc- 
tion reddening vector is appropriate for w 85% of 
PMS stars in a SFR similar to Taurus. Sources for 
which this is not the case (which should be a small 
number), can be pointed out by an unusually larger 
(~30 mag) extinction value. For these sources, the 
extinction reddening vector would rotate in the C- 
M diagram to become nearly vertical, signaling that 
only attenuation of light and no reddening is present. 
In the C-C diagram the vector would decrease in 
magnitude significantly, and even disappear, inval- 
idating the application of the method that we are 
proposing. This includes the most extreme cases of 
a disk viewed edge on, when the light of the star be- 
comes completely extincted and the spectral energy 
distribution is dominated by scattered light from the 
polar regions of the disk. 

2.4. Principal Vectors 

From Figures 1 and 2 we perceive that it is pos- 
sible to represent the IR excess produced by the 
disk (either reprocessing or accretion) as vectorial 
quantity in each C-C and C-M diagram. This vector 
points nearly in the same direction for stars of all 
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masses in the low to intermediate range. For masses 
grater that « 2.5 Mq or spectral types earlier than 
AO the vector starts to rotate, so we regard this as 
the upper limit for the mass. For lower stellar masses, 
the median value of the displacements, represented 
by Z), reproduces the mean effect of the disk. The 
magnitude of the effect is however different for each 
case, but this only implies that an amount ||-D|| of 
disk emission should be interpreted differently for 
each case, depending on the stellar mass. To achieve 
a more realistic approach, we will consider the D 
vectors corresponding to the accretion models that 
were derived in §2.2. Furthermore, we note that the 
direction of D is very different from the reddening 
vector produced by extinction X. This vector seems 
appropriate for star forming regions not affected by 
dense cores of Ay > 25 mag and disks not viewed at 
inclination angles such as Ay > 30 mag that would 
be affected by abnormal extinction. 

Therefore we have a pair of vectors D and X 
for each diagram. We call them Dec and Xcc for 
the C-C D and Den and Xcm for the C-M D. 
The components of the vectors are: I?cc=(0.767, 
0.221), Xcc=(0.6127, 1.064), Dcm={^-01, -1.105) and 
^cm=(1.676, 1.1613) in magnitude units. These two 
sets of vectors are linearly independent so they repre- 
sent an alternative basis for both diagrams. By this 
way we can make new C-C and C-M diagrams us- 
ing D and X as a basis with respect to an arbitrary 
origin. 

So if we have a diskless standard star standard 
without extinction (and at a given distance), it will 
occupy a certain position in each diagram. If we now 
add certain amount for extinction and disk IR ex- 
cess, it will move accordingly in each diagram. On 
the other hand, if we start from a s-f d system, we can 
subtract a certain amount of D and X to recover the 
original location of the star on the diagrams. More- 
over, if we do not know the original stellar mass (or 
spectral type), we can run the same procedure for all 
the possible standard stars (for example, the main 
sequence), and find in which case we get the same 
D and X values for the same standard star on each 
diagram. In the next section, we will show that this 
procedure leads us to identify the star's mass with 
an acceptable error. 

To proceed we need to locate the main sequence 
on the transformed diagrams. And in the case that 
we have young stars we need to place the isochrones 
generated by evolutionary models at a certain age. 
This presents a problem since different models dis- 
agree up to factors of 2 in the masses predicted for 
the same age and spectral type. We will consider sev- 



eral models (see §3.3) but in order to proceed here, 
we will start by using PS99 as a first example of the 
use of the method. To plot the isochrones on the new 
diagrams first we need to know the absolut e JHK 
magnitudes. We use the results of Testi et al. I (|l998h 
for the bolometric corrections made to the PS99 evo- 
lutionary tracks. Later on, we transform the C-C and 
C-M diagrams coordinates to the new basis formed 
by the Principal Vectors (PV). 

3. APPLICATION TO A DYNAMICALLY 
SELECTED SAMPLE OF YSO'S 

To test the PV method, we use a sample of 
Taurus- Auriga (d^l40 pc) PMS with stellar masses 
known by dynamica l methods (TD02: [Schaefer et al 



20031: I Simon et all |2000| ). The majority of the 



Taurus- Auriga objects are intermediate to low mass 
and have a large fraction of binary or multi- 
ple systems. Accord i ng to the study made by 



Kenvon fc HartmannI (|l995r ). the objects are on an 
evolutionary status of no more th an 2 Myr, show ing a 
great amount of IR excesses (see lHaiscbll200l[ ). The 
star forming reg ion is believed to be nearly coeval 
(jHartmannll200ll ). which is one of the advantages in 
studying stellar clusters, that we can assume that all 
the me mbers have an age similar to the age of the 
cluster ( Scalol[l98^ . For this reason first we assumed 
an age of 2 Myr (log(t)=6.3 dex) for all objects. 

Our dynamical sample consists of 14 objects (Ta- 
ble 1) tha t have JHK photoin et ry reported by several 
authors (IStrom et al. _ 1989t Kenvon fc HartmannI 
19951: IWalter et alJIlQSSl ). The extinction reported 
for these objects is generally small {Ay « 30), so 
it is not likely that any of them are viewed edge on, 
and so the PV method can be applied. Some of the 
objects in the sample are known to be multiple sys- 
tems. For this first calculation we took the combined 
photometry of the systems even if they were multiple 
(we did not resolve the components, but see §3.5). 
We applied the PV method to these objects using 
an isochrone of log(t)=6.3 dex. We obtained two di- 
agrams (Figures 3a and 3b) on the new D-X space 
generated by the PV. The axis on both figures are in 
the same scale, so we can compare them directly. All 
the distances on the D-X diagrams are measured in 
extinction units (eU) that corresponds to Ay — 10 
or a disk unity. We can see that if we take an object 
as a reference point on one graph and we put the 
other on top of it, there is a pair of mass values on 
the isochrones that superimpose or are very close to- 
gether. That is, we are finding the mass without disk 
or extinction, and from the vectors, the value of the 
direct contribution of the disk excess and interstellar 
extinction. 
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TABLE 1 

JHK PHOTOMETRY VALUES FOR 
TAURUS-AURIGA OBJECTS THAT HAVE 
MASS VALUES CALCULATED BY 
DYNAMICAL METHODS, IN COLUMN 5 WE 
IDENTIFY IF THE OBJECT IS MULTIPLE. 



object 


J 


H 


K 


Multiple? 


BP Tau^ 


9.30 


8.42 


8.05 


No 


CY Taui 


9.76 


8.90 


8.42 


No 


DL Tau^ 


9.69 


8.77 


8.12 


Yes 


ub iau 


8.32 


7.40 


6.81 


Yes 


DM Tau^ 


10.41 


9.70 


9.45 


No 


FO Tau2 


10.33 


9.34 


8.81 


Yes 


FS Tau2 


10.66 


9.14 


7.74 


Yes 


GG Taui 


8.79 


7.85 


7.25 


Yes 


GM AuP 


9.37 


8.73 


8.48 


No 


LkCalS^ 


9.51 


8.68 


8.22 


No 


UZ Tau El 


8.45 


7.60 


7.02 


Yes 


V397 Aur3 


9.21 


8.46 


8.27 


Yes 


V773 Tau A^ 


7.63 


6.83 


6.48 


Yes 


ZZ Taui 


9.52 


8.78 


8.54 


No 



2 lKenvon fc Hartmaniil lll995l') 
3 lWalter et"aH l|l988h 

We point out that the pair of mass values on the 
two isochrones are separated by a distance, and this 
translates to an uncertainty value of the method. So, 
to calculate the mass of the object, we select the 
mass pair with the "minimum difference" or "error 
measure" (ErrM) between them. 

We can see this more clearly if we plot the ErrM 
vs. mass for each object. In Figure 4 we can ap- 
preciate that the minimum ErrM value corresponds 
to the stellar mass. The mass values are selected 
from a set of given values (0.1, 0.2, 0.4, 0.6, 0.8, 
1, 1.2, 1.5, 2 and 2.5 Mq) that conform the reported 
isochrones by PS99. We are not doing any interpo- 
lation to improve the mass precision. This defines 
the uncertainty on our mass determination. Further- 
more, PS99 did not list temperature and luminosity 
values uniformly with respect to age for each consid- 
ered mass. Because of this, it was necessary to take 
an age dispersion of log{t) = 6.3 ± 0.05 dex, that 
corresponds to an uncertainty value in age of ~ 0.3 
Myr. In the cases where the same mass value ap- 
pears several times on an any given age interval, we 
calculated the ErrM for each one and then obtained 
the mean value. The variations with respect to the 
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Fig. 3. D-X diagrams showing the location of the 14 
Taurus- Auriga objects on the new space made by the 
FY. The numbers represents the masses on a log(t)=6.3 
dex isochrone. The objects are identified by their names. 
The axis units values are extinction units (eU). 

mean values are shown on Figure 4 as vertical bars 
on each point. 

3.1. Comparing the PV method to dynamical mass 
values 

In the last section we obtained the mass values 
for 14 objects applying the PV method. Now, we 
compare these PV masses with the values found by 
dynamical methods, using orbital inform ation of bi- 
nary systems fTD02: ISchaefer et al. II2OO3I) or circum- 
stellar disk rotation kinematics I Simon et al. |[2oo3). 
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GMAur 
mass = 0.6& 

- 2 



UZTauE 
mass=1.5& 

- 2 



PV Mass (Mq) 

Fig. 4. Mpv mass calculated for 14 objects at a given 
age of log(t)=6.3±0.05 from the Taurus-Auriga associa- 
tion. Each graph is labeled by the name of the object and 
calculated mass given in Mq. ErrM represents the min- 
imum difference between the vectors length. The units 
are extinction units (eU). 




1 2 3 

dynamical mass (Mg) 

Fig. 5. Mpv vs Mdyn for the objects on a given age. 

Table 2 (see also figure 5) shows this comparison. 
An overall difference of 29% with dynamical masses 
can be appreciated. However, for most cases the mass 
value fits the dynamical value within the dynamical 
error margins. Of the 14 objects, the uncertainty is 
less than la in 6 cases, less than 1.5cr in 4 more and 
less than 3cr in 2 more. The last 2 cases are V397 



TABLE 2 

RESULTS OBTAINED WITH THE PV METHOD 
AT A GIVEN AGE OF LOG(T)=6.3±0.05 DEX, 
FOR THE TAURUS- AURIGA OBJECTS. Mpv, 
DISK VECTOR MAGNITUDE VALUE, 
INTERSTELLAR EXTINCTION, AND 
MINIMUM ERRM VALUE. 

Object Mdyn ±err Mpyierr D X ErrM 



BP Tau 
CY Tau 
DF Tau 
DL Tau 
DM Tau 
CM Aur 
FO Tau 
FS Tau 
GG Tau 
LkCal5 
UZ Tau E 
V397 Aur 
V773 Tau 
ZZ Tau 



1.242 
0.552 
0.82* 
0.722 
0.552 
0.842 
0.774 
0.78* 
1.54" 
0.972 
1.31 2 
2.26^ 
3.20* 
O.8II 



+0.25 
-0.32 

0.33 
0.24 
0.11 
0.03 
0.05 
0.25 
0.25 
0.78 
0.03 
0.08 
0.21 

0.71 

+0.44 

-0.25 



l.OiO.l 
0.60±0.1 
2.0±0.25 
0.60±0.1 
0.40±0.1 
0.60±0.1 
0.60±0.1 
0.80±0.1 
1.5±0.15 
0.80±0.1 
1.5±0.15 

l.OiO.l 
2.5±0.25 
0.60±0.1 



0.192 
0.112 
0.303 
0.110 
0.067 
0.045 
0.275 
0.569 
0.247 
0.118 
0.141 
0.129 
0.596 
0.038 



0.178 
0.320 
0.431 
0.553 
0.004 
0.136 
0.227 
1.172 
0.452 
0.301 
0.519 
0.023 
0.011 
0.084 



0.045 
0.004 
0.078 
0.053 
0.127 
0.048 
0.126 
0.034 
0.019 
0.080 
0.062 
0.044 
0.200 
0.074 



O lGuilloteau et al.l l ll99gt)IWhite et al.l lll999l') 
l lSchaefer et~all (120031) 
2i Simon et all l|2000t ) 
3 tSteffen et al.l l|200ll) 
4 lTamazian et al. 1 1I2OO2D 

Aur (AM = 1.26) and DF Tau (AM = 1.18). We 
will explore possible causes for the discrepancies. 

There is a controversy on distance on DF Tau. 
Hipparcos measurements ( Favata et alj 1998) place 
it at a distance between 31-52 pc. iLamzin et al. 
(l200Th based on HST measurements debate this 



value and place it at 70 pc. However. ISchaefer et al 



(|2003l ). determine the system's dynamical mass to 
be 0.9 Mq assuming a distance of 140 pc, the same 
distance used by TD02. On the other hand, if the 
distance were as small as 70 pc, DF Tau would not 
be a Taurus member (extension of '^20 pc), and the 
dynamical mass would change as . So we will take 
the fiducial distance of 140 pc for DF Tau. V397 
Aur (NTT 045251+ 3016) is a spectroscopic binary 
( Walter et al.lll988h and is classified as a naked T 



Tauri . The distance is also a problem, iSteffen et al 



(|200ll ) calculated their dynamical masses setting the 
system at 145±8 pc. If we took DF Tau out of the 
comparison we would obtain a 20% agreement with 
dynamical masses. 




Fig. 6. Mpv vs Mdyn for the objects at the age of better 
match (hsted). 



We remark the fact that the calculation has been 
made taking a given age (coeval). However there 
are indications that star clusters have an age dis- 
persion. For the Taurus region several authors have 
mentioned that this varia ti on is ±0.3 dex on ag e 
(|Hartigan &: KenvonI 120031 : IPalla &: Stalheij l2002l ). 
This can be another issue to take into account to try 
to explain the observed differences. To prove this last 
hypothesis, we will apply the PV Method backwards, 
that is, since for these objects we know the dynami- 
cal mass value, we will find the isochrone that gives 
the minimum ErrM for that mass. In the next sec- 
tion we will address the question of obtaining mass 
and age simultaneously. The outcome age values are 
listed in Table 3. (See also Figure 6). There we com- 
pare them to the ages estimated by other authors. 
We appreciate in general a certain dispersion in the 
ages determined by other authors and in particu- 
lar for the age of DF Tau, all seem to agree that 
the age is significantly younger than log(t)=6.3 dex. 
We get a value of log(t)=5.55 dex, which is similar 
to the one previously determined bv lPalla fc Stalheij 
( 20021 ) that assigned it an age of log(t)=5 dex, so 
probably DF Tau is younger than the cluster's age. 
For V397 Aur the age value found was log(t)=6.5 
dex that corresponds to the upper limit for which we 
have data values on the evolutionary track for masses 
more than 2.5 Mq (see also §3.2), so this object is 
out of the boundaries of the restrictions imposed by 
the method (more massive and older). 



Fig. 7. ErrM as a function of mass and age for each 
object in the sample. The solid line represents the region 
where ErrM has a minimum (<0.3 eU with a step of 0.05) 
value and this constitute possible solutions for mass and 
age. The dashed line corresponds to 0.3 eU<ErrM<0.6 
eU with a step of 0.1 and the dotted lines to an ErrM 
value greater than 0.6 eU (step of 0.25). 

3.2. Mass and Age 

We have seen that we can obtain different mass 
values with the PV method if we change the age. 
We want to explore if it is possible to obtain the 
age and mass simultaneously for each object. That 
is to develop a method to determine the two quan- 
tities. However, this solution involves 4 quantities 
(mass, age, size of the disk and extinction redden- 
ing) and we only have 3 photometric measurements 
(J, H and K). So it is expected that the solution 
is not unique, and dependent of one of the parame- 
ters. We apply the PV method for a grid made by 
age values and visualize the results as a contour map 
with the ErrM error as a function of two variables 
(mass and age). We can see how the minimum value 
changes as a function of mass and age. These maps 
are shown in Figure 7. We can detect the presence 
of a region on the top left of each panel (continuous 
solid line) where a canyon like shape is formed. This 
region is generated by all the minimum ErrM that 
corresponds to all the possible acceptable solutions 
for mass and age for the object that are consistent 
with the measurements of J, H and K. 

For DF Tau and V773 Tau (Figure 7) we can bet- 
ter appreciate that the contour figures are slightly 
different from the rest of the objects. The solid line 
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TABLE 3 

AGES DETERMINED BY DIFFERENT AUTHORS, COMPARED WITH PV OBTAINED VALUES. 



Object 


Mdyn 


Mpv 


PS02 


WGOl 


SOO 


StOl 


H95 


TOl 


G04 


Agepy±0.12 


BP Tau 


1.24 


1.0 


6.25 


6.16 


6.3-7.0 




5.78 






6.5 


CY Tau 


0.55 


0.60 


6.25 


6.48 


6.3-6.7 




6.27 






6.25 


DP Tau 


0.82 


1.0 


5.0 








3.74 






5.5 


DL Tau 


0.72 


0.80 


6.4 


<6 


6.30 




5.78 






6.5 


DM Tau 


0.55 


0.60 


6.44 


6.80 


6.70 






6.4 


6.4 


6.75 


FO Tau 


0.77 


0.80 


5.6 














6.5 


FS Tau 


0.78 


0.80 


6.44 














6.25 


GG Tau 


1.54 


1.5 


6.1-6.36 




6.0-6.3 




4.82 




5.8-6.2 


6.25 


GM Aur 


0.84 


0.80 


6.98 


6.88 


6.47 




5.06 


6.25 


6.1-6.25 


6.5 


LkCa 15 


0.97 


0.80 


6.71 


6.80 


6.4-6.7 






7.07 


6.9-7.07 


6.75 


UZ Tau E 


1.31 


1.5 


5.3 








3.60 






6.0 


V397 Aur 


2.26 


2.5 


6.52 






6.78 








6.5 


V773 Tau 


3.20 


3.0 


6.44-6.6 














4.0 


ZZ Tau 


0.81 


0.80 


5.6 














6.5 



PS02: Palla fc Stalher (2002); G04 : Greaves et al. (2004); H95:H artigan et aL (1995) 
WGOl: White fc Ghez (2001); S00:Simon et al. (2000): StOLSteffen et aL (2001) 



zone is located at a lower position, indicating that 
these objects are younger than the rest, with ages 
less than 10^ yr for any given object's mass. On the 
other hand, DM Tau shows the solid line zone lo- 
cated on a higher zone indicating an older age. For 
V397 Aur, the zone extends to an older region as the 
mass is incremented, getting out of the range when 
the mass is larger than ~2M0, restricting us from 
finding the real mass as we already mentioned. 

3.3. Choosing a different evolutionary models 

We are bound to evolutionary track models, 
because those are used as a reference points by 
the PV method in the D-X space. We explore 
the effect on mass determination if we take dif- 
ferent evolutionary track models. In order to ad- 
dress this point, we used several evolutionary mod- 
els such as DM97, BCAH98, SOO and F^qi, each 
one wi th bolometric correct i ons f ro m different au- 
thors ( Kenvon fc HartmannI 1995 ; Leieune et al. 
19981: ISiess et al.lll997h and we calculate the resulting 
Mpv masses using the coeval age of log(t)=6.3 dex. 
Then we compare the results for each model with 
dynamical masses in table 4. In all the cases we 
took DP Tau out of the comparison for the reasons 
explained above. For the mass values obtained with 
DM97 evolutionary models for a mass range from 
0.017 to 2.5 Mq we obtained a difference of 23%. 
For the BCAH98 tracks we did not take into account 



V773 Tau and V397 Aur because they have dynami- 
cal masses greater than the superior mass limit that 
is possible to determine with them (0.07 to 1.4 Mq), 
obtaining a value of 33% of variation with respect to 
dynamical masses. 

The SOO evolutionary models consider an accre- 
tion of Mace — I.O.tIO^^ M^jlyr over a period of 
5a;10^ yr onto the central star that produces an ex- 
tra luminosity and an increase of mass. For masses 
smaller than 0.4 Mq a problem exists because the 
accretion processes directly affects the object, caus- 
ing that an object with an initially mass of 0.2Mo 
will appear on the evolutionary track in a locus of 
a O.5M0 star of the corresponding age. After apply- 
ing the method we obtained a difference of 34% for 
the model with overshooting and 45% without over- 
shooting. If we use the F^Ol for a Z=0.02 and mass 
range from 0.40 to 2.6 Mq we obtained a variation 
of 23% with respect to dynamical. 

We thus see that the evolutionary models pro- 
posed by DM97, F^, and PS99 give all similar val- 
ues that compare to dynamical masses at the ^ 20% 
level. Furthermore, DM97 and compare also rea- 
sonably good 23% ) to PS99 tracks that we em- 
ployed through §3. 

Therefore, at this level of certainty we consider 
any of these models as equivalent, and we will con- 
tinue to use the PS99 tracks. However, we must em- 
phasize that the PV method is very sensitive to the 
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particular evolutionary model employed, and since 
this topic is not yet settled, people should remain 
aware of further developments in this area. 



3.4. Uncertainties 

Several other factors have to be considered to put 
out the method into context, such as the effect of dis- 
tance, age uncertainties. In this section we estimate 
these and other factors and quantify how much the 
PV mass value for an object is modified in each one 
of the listed conditions. We first explore variations 
in the distance to the object. That is to say, if we as- 
sume that the region is not at 140 pc. In order to see 
this, we made an analysis varying the distance in ± 
10 pc (7%), and we applied the PV method with the 
new parameter to the sample objects. The masses 
obtained yield a median difference of 22% with re- 
spect to the original PV value at 140 pc. To examine 
the uncertainty that appears when we assign a dif- 
ferent age, we took PS99 tracks and we varied the 
age in ± 0.1 dex (22%) and obtained a 18% of vari- 
ation in the Mpv determination with respect to the 
coeval age of log(t)=6.3 dex. 

Another source of uncertainties could be due to 
the principal vectors themselves. If a different value 
of R is assumed, the slope of X vector will only 
change in the C-M D. But it will remain the same 
in the C-C D. To quantify the magnitude of the PV 
mass produced by this change of slope we tried val- 
ues of R (3.2 and 4) instead of the normal value of 
R=3.1 for the ISM. The resulting PV masses differ in 
1.5% and 37% respectively. Therefore, small change 
in the value of R has no significant effect on mass 
but a larger change such as required for dense cores 
(R~4) should be carefully considered. 

Next, we want to evaluate the change in PV mass 
when the slope of the excess vector D varies. We 
use th e slope value of 0.58 obtained bv iMever et al 



(|l997f ) for the locus of the T Tauris as the slope value 
for the excess vector D. We assigne d the magnitude 



|D|| obtained through the models of iD'Alessio et al. 



(j2005D of \\Dcc\\ = 0.798. With these, values we ap- 
plied the coordinate transformation to the D-X space 
at a given age of log(t)=6.3 dex. We needed two dia- 
grams to find the mass values but lMever et al.l (|l997( ) 
obtained his T Tauris locus only in C-C D. So, we 
decided to use the C-M D slo pe for the excess vector 
that was obtained with the iD'Alessio et al. 



( 2n05[ ) 



accretion models, which has a slope of -1.09. When 
we changed these parameters, we obtained a mass 
variation of 12% with respect to the values obtained 
in S3.1. 



TABLE 4 

MASS OBTAINED USING SEVERAL 
EVOLUTIONARY MODELS. AGE 
LOG(T)=6.3±0.05 DEX. THE UNITS ARE IN 

SOLAR MASSES. aovN IS THE MEAN 
RELATIVE DIFFERENCE TO DYNAMICAL 
MASS. 



^Djeci 










A-IDM97 




iVlYY 


BP Tau 


1.24 


1.0 


1.0 


1.0 


1.30 


0.80 


1.10 


CY Tau 


0.55 


0.6 


0.60 


0.6 


0.80 


0.45 


0.60 


DF Tau 


0.82* 


2.0* 


1.89* 


3.5* 


2.5* 


1.1* 


5.0* 


DL Tau 


0.72 


0.6 


0.60 


0.6 


0.80 


0.5 


0.60 


DM Tau 


0.55 


0.4 


0.13 


0.1 


0.40 


0.30 


0.5 


FO Tau 


0.77 


0.6 


0.16 


0.5 


0.60 


0.40 


0.5 


FS Tau 


0.78 


0.8 


0.25 


0.2 


1.0 


0.57 


0.80 


GG Tau 


1.54 


1.5 


1.29 


1.3 


2.0 


1.15 


1.80 


GM Aur 


0.84 


0.6 


0.5 


0.5 


0.80 


0.45 


0.40 


LkCa 15 


0.97 


0.8 


0.25 


0.6 


0.80 


0.5 


0.80 


UZ Tau E 


1.31 


1.5 


1.39 


1.4 


2.0 


1.40 


1.90 


V397 Aur 


2.26 


1.0 


0.80 


0.8 


1.2 


0.60** 


1.0 


V773 Tau 


3.2 


2.5 


6.0 


2.5 


3.0 


1.10** 


2.90 


ZZ Tau 


0.81 


0.6 


0.60 


0.6 


0.80 


0.5 


0.60 



(Tdyn (%) 20.29 44.9 34.6 22.94 33.64 23.6 

l:without overshooting 
2:overshooting 

*:Not taken into account in the cr^yn 
**:Out of mass range defined by the ' 



value, see text. 

BCAH98 model (0.07 to 1.4 Mq). 



3.5. Binaries 

Many of the objects in our list belong to multi- 
ple systems and we have analyzed them as a single 
object. So the next question emerges. Can we apply 
the method to the components of a system and find 
the individual masses in a consistent way? To an- 
swer this, we need to know individual photometry for 
the components and apply the PV method for each 
one an determine the mass for a given age. However, 
we have dynamical mass values only for V773 Tau, 
the other systems h ave their masses calc ulated by 
evolutionary tracks ( White fc Ghejl200ll ) (hereafter 
WGOl), and are taken so that the sum of the compo- 
nents match the mass values obtained by dynamical 
methods. We decided to apply the PV method to 
the systems components and compare these values 
with mass values determined by evolutionary tracks. 
Only GG Tau, FO Tau, FS Tau, and V773 Tau, meet 
the requirements as they have complete component 
photometry and determined masses (see table 5). 
For the preliminary determination we fixed the age 
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TABLE 5 

MULTIPLE SYSTEMS PHOTOMETRY. 



Components 



AJ H AH K AK 



GO Tau Aai 
GGTau Abi 
GG Tau Ba^ 
GG Tau Bbi 
FS Tau A2 
FS Tau B2 
FO Tau A2 
FO Tau B^ 
V773 Tau AB^ 
V773 Tau C^ 



9.24 0.09 8.27 0.09 

10.12 0.02 9.07 0.08 
11.48 0.16 10.63 0.15 
13.16 0.12 12.38 0.06 
10.85 0.0 9.32 0.0 
12.66 0.0 11.17 0.0 

10.13 0.0 9.35 0.0 
10.95 0.0 9.59 0.0 
7.77 0.05 7.03 0.03 
10.13 0.09 8.91 0.15 



7.73 0.05 
8.53 0.08 
10.20 0.12 
12.01 0.22 
7.78 0.0 
10.03 0.0 

8.76 0.0 
9.14 0.0 

6.77 0.09 
8.09 0.31 



1 : White et al. (1999); 2 :Woitas et al. (2001) 

to log(t) =6.3 dex for al l the s ystems, following the 
work of iHartigan et al" (1994) where they studied 
39 binary systems in Taurus-Auriga and Orion re- 
gion, finding that 2/3 of the sample we re coeval. GG 



Tau i s a quadruple interactive system (jLeinert et al 



19911 ) of two binary pairs identified as GG Tau A and 
GG Tau B separated by a distance of 10." 1 (1,414 
AU). For the GG Tau A components the separation 
is 0."25 (35 AU)( named GG Tau Aa and GG Tau 
Ab), (see mass values in table 6). The GG Tau B 
components are separated by l."48 (207 AU). In Ta- 
ble 4_we_sh£w_the_JFIK ph otometric va l ues re ported 
bv lWhite et all ([l99i) and lGhez et aD (| 19971) . 

The FS Tau system (Haro 6-5 A) has two ob- 
jects separated by a distance of 0."25, and the sec- 
ondary component isj^ times weaker than the pri- 
mary (jChen et al.lll990l ). TD02 determined for the 
first time the dynamical value for the system's mass. 

V773 Tau is a multiple s y stem, a t dist ance of 
148 pc (|Lestrade et al.lll999l ). IWeltvl (|l995t ) deter- 
mined it as triple system with a very close pair 
(0.3AU) and the other component separated at 0." 15 
(^22 AU). But in 2 003, several authors r eported 
a fou r th component (iDuchene et al. 120031: IWoitad 
l2003h . iDuchene et al.l ( 2003 ) renamed the closer bi- 
nary pair (System V773 Tau AB) as V773 Tau A 
and V773 Tau B and also V773 Tau C for the one 
located at 0."15 and V773 Tau D for the new com- 
ponent that was cataloged as an IR Source. 

In Table 6 and Figure 8 we show the resulting 
mass values with the PV Method at a given age 
of log(t)=6.3 dex. In general we can see that the 
method converges for most of the sample objects. 
Only for the case of GG Tau Bb we can not find 
the minimum because we are mass limited. We pre- 



C 



o 



1 Z 30 1 3 30 1 3 30 1 3 



1 3 3 



PV Mass(Mg) 

Fig. 8. Determined PV mass for the components of GG 
Tau, FS Tau , FO Tau and V773 Tau. The PV Mass val- 
ues are shown on the top right margin under the objects 
name on each graph. 

diet for this component a mass of <0.1 M© having 
an ErrM of 0.84 eU. For the 10 of the objects listed 
in table 6 (Figure 9), 7 coincide better than 1.5cr 
according to track mass uncertainties and 1 more 
is closer than 2ct. Only 2 cases are >3cr away values 
(GG Tau Aa and FO Tau B). We ran the PV method 
with a varying age to try to find if we can explain 
these discrepancies by finding the objects age. We 
show the results in table 7 and Figure 10. We can 
see that we can obtain a fitting value within Ict, if 
the ages of GG Tau Aa and FO Tau B arc log(t)=6.0 
dex. The age determined for GG Tau Aa matches the 
estimated by WGOl and for FO Tau B we are within 
the uncertainties. 

4. CONCLUSION 

We explored the effects of excess emission from 
disks on the JHK-NIR colors on C-C and C-M di- 
agrams. After analizing reprocessing and accretion 
disk models, a pattern in J-H vs. H-K and K vs. 
J-K diagrams was found that gave origin to a set of 
principal vectors (PV), originated by the IR excess 
produced by a disk and the interstellar extinction. 
We applied a space transformation to verify the pos- 
sibility of determining the mass of an object at a 
given age. 

To test our PV method we applied it to a sam- 
ple of 14 Taurus-Auriga PMS objects, first without 
resolving the components if they were multiple, with 
the purpose to determining the object's mass. We 
found a fair agreement for 12 cases when we used a 
coeval asumption. For the other 2 cases, in particu- 
lar for DF Tau, we showed that an agreement with 
the dynamical mass is reached for a different age, as 
it had been already pointed out by other authors. 
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1-->GGTauAa _ 
>QGTauAb 
>QGTauBa 
GGTauBb 
FOTauA 
FOTauB 
FSTauA 
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>V773TauAB 
10-->V773TauC 



1 2 

track mass (Mg) 



> 



GGTauAa-at-6dex 
GGTauAb-at-6.25dex 
GGTauBa-at-5.75dex 
GGTauBb-at-6.75dex 
FOTauA-at-6.75dex 
F0TauB-at-6dex 
FSTauA-at-6.25dex 
FSTauB-at-6.25dex 
V773TauAB-at-5.5dex 
10-->V773TauC-at-6dex 




track mass (Mg) 



Fig. 9. Mpv vs Mxrack for the system components 
of GG Tau, FO Tau and FS Tau. The masses for the 
sys tem V773 Tau are d ynamical determinations made 
by iDuchene et all Jixil ). We assume a coeval age of 
log(t)=6.3±0.05 dex for all the systems. 

TABLE 6 

MASS DETERMINATION FOR MULTIPLE 
SYSTEMS COMPARED TO GIVEN MASS BY 
EVOLUTIONARY TRACKS (WGOl). PV MASS 
VALUES SHOWN ON COLUMN 3. COLUMNS 4, 
5 AND 6 DISPLAY THE CALCULATED DISK 
AND INTERSTELLAR EXTINCTION 
VECTORS AND ERRM. 



Object 


Mrrocfcierr Mpyierr 




D 


X 


ErrM 


GG TauAa 


0.76±0.09 


1.0 ±0.1 





194 


0.294 





150 


GG TauAb 


0.68±0.02 


0.60±0.1 





257 


0.172 





174 


GG TauBa 


0.12±0.03 


0.20±0.05 





148 0.0191 


255 


GG TauBb 


0.042±0.019 


0.10±0.05 





285 


0.462 





844 


FO TauA 


0.60±0.09 


0.40±0.1 





004 


0.378 





296 


FO TauB 


0.59±0.09 


O.SOiO.l 





648 


0.236 





188 


FS TauA 


0.53±0.08 


0.60±0.1 





497 


1.240 





270 


FS TauB 


0.20±0.08 


0.20±0.5 





599 


0.556 





217 


V773 TauAB 


2.6i±0.2 


2.50±0.25 





590 


0.266 





175 


V773 TauC 


0.7i±0.1 


O.SOiO.l 





374 


0.466 





187 



1: dynamical mass by Duchene et al. (2003) 



We also applied the method to 4 multiple sys- 
tems and determined the mass of their individual 
components, getting a better than 1.5ct discrepancy 
for most of the components when we compared with 
evolutionary track mass determinations. 



Fig. 10. Mpv vs Mtrack to determine ages for the sys- 
tems components at given mass that fits the mass value 
determined by evolutionary tracks. The top left labels 
displays the name of an object and the respective age 
determinated by the PV method. 

TABLE 7 

AGES (AGEpv) DETERMINED BY THE PV 
METHOD FOR A GIVEN EVOLUTIONARY 

TRACK MASS (MTi?Ac A') FOR THE 
COMPONENTS OF THE FOUR MULTIPLE 
SYSTEMS. 

Object Af^racfe Age^idex) Agepv±0.l2 



GG Tau Aa 


0.76 


0.80 


6.0±0.35 


6.0 


GG Tau Ab 


0.68 


0.60 


6.03±0.25 


6.25 


GG Tau Ba 


0.12 


0.10 


6.26±0.39 


5.75 


GG Tau Bb 


0.042 


0.10 


<6.0±0.3 


6.75 


FO Tau A 


0.60 


0.40 


6.27±0.25 


6.75 


FO Tau B 


0.59 


0.40 


6.30±0.25 


6.0 


FS Tau A 


0.53 


0.80 


7.89±0.3 


6.25 


FS Tau B 


0.20 


0.20 


7.23±0.38 


6.25 


V773 Tau AB 


2.6 


2.5 


6.74±0.17 


5.5 


V773 Tau C 


0.7 


0.6 


<6.0±0.3 


6.0 


l:BarafFe et al. 


(1998) 









2 : White fc Ghez (20011 

Thus, we have shown that this method works for 
a particular kind of PMS objects, that is class II 
classical T Tauri stars. The method still has to be 
tested on youger (class I) or older systems. Currently 
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we can apply the PV method only to a mass range of 
0.1 to 2.5M0. Evolutionary tracks for a wider range 
of masses would be useful to expand the lower mass 
limit of the method. For the upper mass limit, we 
have the additional problem that the disk vector 
shows a steeper slope in the C-C D. Accretion models 
will be needed to further investigate if an apropiate 
set of vectors for higher stellar masses can be found. 
Also, the method needs to be tested in regions with 
abnormal extinction Ay > 25 or R higher than nor- 
mal, where an appropriate extinction vector would 
has to be calculated. For the small fraction of stars 
that would be viewed nearly edge on, gray extinction 
will appear and the method will fail. However these 
cases should be sorted out by their unusually larger 
extinction. 

Variations in the calculated mass for distance and 
age uncertainties have been also evaluated. We con- 
sidered reasonable variations in distance (7%) and 
age (22%) to the objects in our Taurus sample and 
found a variation in the Mpv masss of 20%. This 
value is similar to the 20% residual error of the PV 
method with respect to dynamical masses, sugesting 
that most differences can be explained as an age and 
distance dispersion effect. 

Although we mainly used PS99 tracks, we also 
considered different evolutionary models. The results 
differed from 20% to 45% with most around 22% 
with respect to dynamical masses. The worst com- 
parision is with SOO overshoting models, that have 
the effect of significant accreted mass on some ob- 
jects. The better agreement is with PS99 models. 

As a result of this first analysis, we conclude that 
the PV method is a promising tool to obtain the 
masses of low to intermediate PMS stars from JHK 
photometry. 

For a future work, the method needs to be tested 
with larger sets of data. Either with known dynami- 
cal masses, althougth the advantage of same age and 
distance would be compromised, or exploring its per- 
formance with other clusters with known IMF. 

T.A.L.Ch. was partially supported by CONA- 
CyT project No.36574-E. 

REFERENCES 

Baraffe, I., Chabrier, G., AUard, F., & Hauschildt, P. H., 

1998, (BCAH98), A&A, 337, 403 
Baraffe, I., Chabrier, G., AUard, F., & Hauchildt, P.H., 

2002, A&A, 382, 563 
Bessell, M.S. & Brett, J.M., 1988, PASP, 100, 1134 
Brandl, B.; Brandner, W.; Eisenhauer, F.; Moffat, A. F. 

J.; Palla, F.; Zinnecker, H., 1999, A&A, 352, 69. 



Briceno, C.; Luhman, K. L.; Hartmann, L.; Stauffer, J.R.; 

Kirkpatrick, J.D.; 2002, ApJ, 580, 317. 
Bontemps, S., Andre, P., Keias, A. A., Nordh, L., Olofs- 

son, G., Huldtgren, M., Abergel, A., Blommaert, J., 

Boulangcr, F., Burgdorf, M., and 12 coauthors, 2001, 

A&A, 372, 173. 
Calvet, N., Patino, A., Gladis Magris, C. & D'Alessio, 

P., 1991, ApJ, 380, 617 
Calvet, N., Gladis Magris, C, Patino, A. & D'Alessio, 

P., 1992, RevMoxAA, 24, 27 
Chen, W.P., Simon, M., Longmore, A.J., Howell, R.R. & 

Benson, J.A. 1990, ApJ, 357, 224 
Chiang, E.I. & Goldreich, P. 1997, ApJ, 490, 368 
Chiang, E.I. & Goldreich, P. 1999, ApJ, 519, 279 
Cohen, M. & Kuhi, L., 1979, ApJS, 41, 743 
D'Antona, F. & MazziteUi, I. 1994,ApJS, 90, 467 
D'Antona, F. & MazziteUi, I. 1997, (DM97), 

Mem. Soc. Astron. Italiana, 68, 807 
D'Alessio, P., Canto, J., Calvet, N., & Lizano, S., 1998, 

ApJ, 500, 411 

D'Alessio, P., Mcn'n, B., Calvet, N., Hartmann, L., & 
Montesinos, B., 2005, RevMexAA, 41, 61 

D'Alessio, P., Calvet, N., Hartmann, L., Lizano, S., & 
Canto, J., 1999, ApJ, 527, 893 

D'Alessio, P., Calvet, N., & Hartmann, L., 2001, ApJ, 
553,1,321 

Dahm, S. E., 2005, AJ, 130, 1805. 

Duchene, G., Chez, A.M., McCabe, C, & Weinberger, 

A.J., 2003, ApJ, 592, 288 
DuUemond, C. P.; Dominik, C; Natta, A.,2001, ApJ, 560, 

2,957 

Favata F., Miccla, G., Sciortino, S., & D'Antona, F., 

1998, A&A, 335, 218 
Figuoredo, E., Blum, R. D., Damineli, A., & Conti, P. S., 

2002,AJ, 124, 2739. 
Ghez, A.M., White, R.J., & Simon, M., 1997, ApJ, 490, 

353 

Greaves, J.S. 2004, MNRAS, 351, L99 
Guilloteau, S., Dutrey, A., & Simon, M., 1999, A&A, 348, 
570 

Johns-KuU C, & Valenti, J., 2001, ApJ, 561, 1060 
Haisch, K.E. Jr., Lada, E.A., Lada, C.J. 2001, ApJ, 553, 
153 

Hartigan, P., Strom, K., & Strom, S., 1994, ApJ, 427,961 
Hartigan, P., Edwards, S.,& Ghandour, L., 1995, ApJ, 
452, 961, 736 

Hartigan, P., Kenyon,S., Hartmann, L., Strom, S., Ed- 
wards,S., Welty, A. & Stauffer, J., 1991, ApJ, 382, 
617 

Hartigan, P. & Kenyon,S., 2003, ApJ, 583, 334 
Hartmann L., 2001, ApJ, 121, 1030 
Herbig, G. H., & Dahm, S. E., 2006, AJ, 131, 1530. 
Hillenbrand, L.A., Strom, S., Vrba, F.J. & Keene, J., 

1992, ApJ, 397, 613 
Hillenbrand, L.A., Carpenter J.M., 2000, ApJ, 540, 236 
Hillenbrand, L.A. & White R.J., 2004, ApJ, 604, 741 
Isella, A. & Natta, A., 2005, a, 438, 3, 899 
Kenyon, S. & Hartmann, L., 1987, ApJ, 323, 714 



16 



LOPEZ-CHICO & SALAS 



Kenyon, S. & Hartmann, L., 1995, ApJS, 101, 171 

Lada, C.J. & Adams, F., 1992, ApJ, 393, 278 

Lamzin, S.A., Melnikov, S.Y, Grankin, K.N. & Ezhkova, 

O.V., 2001 , A&A, 372, 922 
Leinert, Ch., Haas, M., Richichi, A., Zinncckcr, H. & 

Mundt, R., 1991, A&A, 250, 407 
Leistra, A., Cotera, A. S., Liebert, J., 2006, AJ131, 2571 
Lejeune, T., Cuisinier, F., Buser, R., 1998, 

A&Asupplements, 130, 65. 
Lestrade,J.F., Preston, R.A., Jones, D.L., Phillips, 

R.B., Rogers, A.E.E., Titus, M.A., Rioja, M.J. & 

Gabuzda,D.C. ,1999, A&A, 344, 1014. 
Levine, J.L., Steinhauer, A., Elston, R.J., Lada, E.A., 

2006, AstroPH 4315. 
Liu, Wilson M., Meyer, Michael R., Cotera, Angela S., 

Young, Erick T., 2003, AJ, 126, 1665. 
Luhman, K.L., Rieke, G.H., Lada, C.J., Lada, E.A., 1998, 

ApJ, 508, 347. 
Luhman, K. L.,1999, ApJ, 525, 466. 
Luhman, K.L., Rieke, G.H., 1999 ApJ, 525, 440. 
Luhman, K.L., 2000, ApJ, 544, 1044. 
Luhman, K.L., Rieke, G.H., Young, E.T., Cotera, Angela 

S., Chen, H., Rieke, M.J., Schneider, G., Thompson, 

R.L, 2000, ApJ, 540, 1016. 
Luhman, K.L., Stauffer, John, R., Muench, A. A., Rieke, 

G.H., Lada, E.A., Bouvier, J., Lada, C.J., 2003, ApJ, 

593, 1093. 

Luhman, K. L., Lada, E.A., Muench, A. A.; Elston, R.J., 

2005, ApJ, 618, 810. 
McCabe, C, Ghez, A.M., Prato, L.,Duchene, G., Fisher, 

R. S., Telesco, C, 2006, ApJ, 636, 932. 
Meyer, R.M., Calvet, N., Hillenbrand. L.A., 1997, AJ, 

114,(1) 288 

Moore, T. J. T., Lumsden, S. L., Ridge, N. A., Puxley, 

P. J., 2005, MNRAS, 359, 2, 589. 
Muench, A., Lada, E.A., Lada, C.J. & Alves, J., 2002, 

ApJ, 573, 366. 
MuzeroUe, J., Calvet, N., Hartmaim, L., D'Alessio, P., 

2003, ApJ, 597,2, L149-L152. 
Natta, A., Meyer, M. & Beckwith S., 2000, (NOO), ApJ, 

534, 838 

Nielbock, M., Chini, R., MUer, S. A. H., 2003, A&A, 408, 
245. 

Oasa, Y., Tamura, M., Sugitani, K., 1999, ApJ, 526, 336. 
Palla, F. 1993, IV-CIWSA, p.43 

Palla, F. & Stalher, S.W., 1999, (PS99), ApJ, 525, 772 
Palla, F. & Stalher, S.W., 2002, ApJ, 581, 1194 
Porras, A., Cruz-Gonzalez, I. & Salas, L., 2000, A&A, 
361, 660. 

Rieke, G.H. & Lebofsky, M.J., 1985, ApJ, 288, 618 
Shakura, N. I. & Sunyaev, R. A., 1973, A&A, 24, 337 
Scalo J.M., 1986, Fund. Cosmic Phys., 11, 1 
Schaefer, G.H., Simon, M., Nelan, E., & Holfeltz, S.T., 

2003, ApJ, 126, 1971 
Siess, L., Forestini, M. & Dougados, C, 1997, A&A, 324, 

556 

Siess, L., Dufour, E. & Forestini, M., 2000, A&A, 358, 
593 



Simon, M., Dutrey, A. & Guilloteau,S., 2000, ApJ, 545, 
1034 

Steffen, A.T., Mathieu, R.D., Lattanzi, M.G., Latham, 
D.W., Mazeh, T., Prato, L., Zinnecker, H. & Loreg- 

gia, D., 2001, AJ, 122, 997 
Strom, K., Strom, S., Edwards, S., Cabrit, S. & Skrut- 

skie, M., 1989, AJ, 97(5), 1451 
Schmidt-Kaler, T. H., 1982, Physical parameters of stars, 

Landolt-Bornstein, New series, Vol VI/2b, Astronomy 

& Astrophysics, Stars and Star clusters, edited by K. 

Schaifer & H.H. Voigt (springer. New York) p. 451 
Tamazian, V.S., Docobo, J. A., White, R.J. & Woitas, J., 

2002, (TD02), ApJ, 578, 925 
Testi L., Palla, F., & Natta, A., 1998, A&A, 133, 81 
Walter, F.M., Brown, A., Mathieu, R.D., Myers, P.C. & 

Vrba, F.J., 1988, AJ, 96, 297 
White, R.J., Ghez, A.M., Reid, I.N., & Schultz, G., 1999, 

ApJ, 520, 811 

White, R.J. & Ghez, A.M., 2001, (WGOl), ApJ, 556, 265 
Welty, A. D., 1995, AJ, 110, 776 

Woitas, J., Leinert, Ch. & Kohler, R., 2001, A&A, 376, 
982 

Woitas, J., 2003, A&A, 406, 685 

Yi, S., Demarque, P., Kim, Y.C., Lee, Y.W., Ree, C.H., 
Lejeune,. T. & Barnes, S., 2001, (Y^Ol), ApJS, 136, 
417 



MASS DETERMINATION FOR T TAURI STARS FROM JHK PHOTOMETRY 



17 



T. Lopcz-Chico and L. Salas: Instituto dc Astrononifa, Univcrsidad Nacional Aiitonoma dc Mexico, Apdo. 
postal 877, 22830 Ensenada, B.C., Mexico (tania@astrosen.unam.nix), (salas@astrosen.unam.mx) 



